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Abstract. We present a mathematical model: dynamical systems over finite 
sets (DSF), and we show that Boolean and discrete genetic models are special 
cases of DFS. In this paper, we prove that a function defined over finite sets 
with different number of elements can be represented as a polynomial function 
over a finite field. Given the data of a function defined over different finite sets, 
we describe an algorithm to obtain all the polynomial functions associated to 
this data. As a consequence, all the functions defined in a regulatory network 
can be represented as a polynomial function in one variable or in several vari- 
ables over a finite field. We apply these results to study the reverse engineering 
problem. 



1. Introduction 

In this paper we introduce the definition of dynamical systems over different finite 
sets (DSF) and we develop its applications to regulatory networks and the Reverse 
Engineering Problem. We consider variables over sets with different numbers of 
elements and we change that to variables over a finite field. 

The justification for considering dynamical systems over different finite sets is 
related with the method Generalized Logical Networks developed by Thomas and 
colleagues, j^lEl^SElEl- The generalized logical networks has a mean consid- 
eration: a variable can have more that two possibilities but always the number of 
possibilities is finite. In addition, the network is described by a function which acts 
over several variables and for each variable there are different number of values. 
These considerations are very important for biologists because it is known that in a 
regulatory network all the variables do not have the same number of states. Here, 
we prove that all of these functions can be considered over a finite field and as a 
consequence of that we can represent them by polynomial functions. In section 2 
we present an algorithm which changes a function over different set of values to a 
function over a finite field. 

In Section 2 we introduce the method to construct functions over a finite field 
using functions defined over finite sets with different number of elements. In Section 
3, we apply the partially defined functions to Reverse Engineering Problem. In 
Section 5, we introduce the definition of Dynamical Systems over different finite 
sets. 
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2. Partially defined functions 

Now, in this section we introduce the mathematical background which will permit 
the application of modelling methods such as generalized logical networks. 

Let Xj — {0, 1, ... ,j — 1} and let Zp be the set of integers modulo p, with p a 
prime number. Suppose that p > j, and we consider a canonical map from Xj to 
the field Zp given by a — > a( mod p). In the following we consider Xj C Zp. Let 
X — [xi, . . . ,Xn) G Zp". We denote the polynomial ring in n variables over Zp by 
Zp[a;i, . . . Xn]. We begin with some definitions. Let D C Zp. 

Definition 2.1. Let S C Zp". Let f : S ^ D be a function. We will call f a 
partially defined function over Zp. 

Example 2.2. Now, let g : X2 x Z3 ^ X2, given by the following table: 



g{xi,x2) 


1 2 




1 


1 1 
1 



Then the table of values of the partially defined function g is the following: 



9ixi,X2) 





1 


2 





1 





1 


1 








1 


2 


* 


* 


* 



X • • • X X„i,^ C Zp", and let D C 'Zp. Since a partially defined 
function f : S D is not a function from Zp" to Zp, we are interested in solving 
the following problem: 

[DF(Zp):] Let S C Zp" and let / : S* be a function. We want a polynomial 
function P : Zp" ^ Zp such that P{x) = f{x), for aU x G 5 C Zp". 

A function P associated to / will be called a polynomial function for f. Now, 
we prove that the problem DF(Zp) can have more than one solution. 

Proposition 2.3. For each function f : S ^ D there is a polynomial P{x) G 
7jp[xi, . . . , Xn], such that f{x) — P{x) for all x G S. The polynomial P can be 
chosen with degree less than or equal to n{p ~ 1) but in general, it is not unique. 

Proof. If /c is a finite field and / : fc" ^ fc is a function then there exists a polynomial 
P in the variables a;i, . . . , Xn, with coefficient in fc, such that f(x) — P(xi, . . . , Xn) 
for all (xi, . . . ,Xn) S fc", [7]. But, in our case we do not have a function from 
/c" — > k, so we will prove that the polynomial function exists associated to the 
partially defined function but it is not unique. 

We will show the idea using the example 12. 21 In the table of / we can complete 
the table in some way. Then there exists a unique polynomial for this table. But 
we can complete the tabic in many ways, so the polynomial function exists but, it 
is not unique. □ 

As a consequence of Proposition 12.31 we have an algorithm which solves the 
problem DF(Zp). Let / : ^ Zp be a function. Let m = \S\ be the cardinality of 
S (m — X]"=i '^i when / : X^i x • ■ • x Xm„ Zp). Now, we write a polynomial P 
in n variables xi, . . . , Xn. P has degree less than or equal to p — 1 in each variable, 
so has degree less than or equal to n{p — 1). We denote P in the following form: 
P{xi, . . . ,Xn) = X^asz " ^aX°' , whcrc a — (ai, . . . , q;„), x°' — x"^ . . . x"" . Now, we 
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evaluate P for all a G 5 and we obtain a system of m linear equations in the 
unknowns ha which always has solutions. The system is the following: 

(/) (a)"6„ = /(a) for all aeS. 

Solving the system using elementary row operations, we finally obtain all the so- 
lutions. In it is proved that the rank of this system is m. Then, there are 
hp-^ , • . . , 6/3„ coefficients of the polynomial P whose are determined in term of the 
free coefficients denoted by 6-^^ , . . . , 6-yp„_,^^ . Now, let Zp^^~^^ [xi, . . . , x„] be the 
subspace of Zp[a;i, . . . , x„] of all polynomials with maximum degree p — 1 in each 
variable and coefficients in Zp. So, we have the following theorem. 

Theorem 2.4. All the polynomial solutions with degree < n(p — I) of the problem 
DF{Zip) are given by a particular solution /o(a?) of (/) plus the subspace 

C/ = {.g e Zp(P-i) [xi,..., x^] \g{a) = 0, Va G S} 
of dimension p^ ~ m. 

Proof. We know by linear algebra that all the solutions of (/) are given by 

fo{x) + b^^gi{x) H h b^^„_^^^gpr,_„,{x) 

where by-^^ , ■ ■ ■ , ^pn_„, G Zp and gi, . . . , gpn^„i G U . Let hi and /i2 be two polyno- 
mial solutions of (/). Then /ii — /i2 G U , so the theorem holds. □ 

3. Reverse Engineering Problem 
over finite sets 

Now, we connect the problem DF(Zp) with the Reverse Engineering Problem 
over Zp. The problem for partially defined functions is equivalent to the following. 

[P(Zp):] Given ai, . . . ,a^ G Zp", b = (fei, . . . G Zp"*, with m < p". Find a 
polynomial P G Zp[xi, . . . , Xn] such that P{sij) — bj for j = 1, . . . , to. 

The problems P{'^p) and DF(Zp) are equivalent. In fact, we only need to take 
S = {ai, . . . ,a„i} and b = (/(ai), . . . ,/(a„)). 

The problem P{Zp) was solved by E. Green in |2. He called the problem P^Zp) 
for (ai, . . . , a™; b) and he proved that if P^Zp) has solutions then the Reverse 
Engineering Problem over Zp has solutions. 

Now, we define the Reverse Engineering Problem over sets with different number 
of elements. Let {kj} be a family of n finite sets where \kj\ = rrij. We denote 
by fc = fci X • ■ • X A:„. Let ri, Ym+i G k. We assume that the vectors 

Vj = {rji, . . . ,rjn) are obtained by experiments (like microarray) and we assume 
that Tj determines Yj+i- Then the Reverse Engineering Problem over k is to find 
a function F — (/i, . . . , /„) : k ^ k such that F{rj) = rj+i for j = 1, . . . , m. But, 
we rewrite: 

[(REP)] The Reverse Engineering Problem over k is to find polynomial functions 
fs'.k—^ks such that 

fsirj) = T-sj+i for j = 1, . . . , TO and s = 1, . . . , n 

Now, we prove that if we can solve DF(Zp), we can solve (REP) and use the 
same algorithm. In fact, if Zp is the field such that p > ruj for all j, we take the 
partially defined functions fs over Zp considering k^ C Zp. Let r = (ri, . . . , r„) G k. 
So, we consider S C Zp" and fs(jj) — rgj+i for j — 1, . . . ,m and s — \, . . . ,n. 
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We have proved the foUowing proposition. 



Proposition 3.1. The reverse engineering problem over set with different number 
of elements has polynomial solutions by Provosition \2.S[ and Theorem[ 

Definition 3.2. 

REP. 



The matrix A — {rj)mxn will be called the matrix of the problem 



Example 3.3. Suppose we have the following data: ri = (1,2,0), r2 = (2,2,1), 
r3 = (1,0,1), r4 ~ (0,1,1), and = (1,1,0). And, we have the additional infor- 
mation: 

(a) the variables {x,y,z} are defined over different finite sets, but we take finite 
fields: cc, y £ Z3 = {0, 1, 2} and z ^ 2,2 — {0, 1}. 

(b) the variable x depends of x and z, y depends of x and y, and z depends of y and 
z. 

The matrix A of the problem is the following: 

/ 1 2 \ 

2 2 1 
A= 10 1 

1 1 
V 1 1 / 

We want polynomial functions /i, /2, and /s, such that F = (/i, /2, /s) and F{Tj) = 
Tj+i for j = 1,2,3,4. 

The additional information (b) means that the functions that we are looking for 
are as follows: 



fi{x,z) 
f2{x,y) 

Using the data we have the table of /i 



aix + a2Z + a^xz + a^x + 052 
+agx^z + O'lxz'^ + flgs^z^ 
60 + bix + b2V + b^xy + b^x'^ + b^y'^ 

+bQX^y + b^xy^ + bgx^y"^ 
Co + ciz + C2y + c^yz + C4y^ + C5Z^ 
+C6y^z + cryz^ + c^y'^z^ 



fi{x,z) 





1 


2 





* 


1 




1 


2 







2 


* 


1 


* 



Using the above table and the algorithm for problem DF(Zp), we obtain a system 
of 4 linear equation with 9 unknown. The matrix of the system of linear equation 
is: 

/ 1 1 1 
12 12 1112 1 
111111111 
V 1 1 1 
Using elementary row operations, we have the following: 

1 + 2a3 + 2a5 + 2aQ + 2a7 + 2ag, — 1 



A, 



2\ 
1 


1 / 



a2 



2a6 + 2as, 



ai ~ 1 + 2a3 + 2a7, ao = as + aj + ae + as- 



APPLICATION OF FINITE FIELDS 



5 



A particular solution of the system is fi — x + z + x'^ . And all the solutions are 
given by /i + 03^1 + 05512 + aegs + arg^ + aggs, where 

gi ^ 1 + 2x + 2z + xz, g2 = 2z + z^, 53 = 1 + 2z + 2x'^ + x^ z 

g4 = l + 2x + 2z + xz^, and 35 = 1 + 2z + 2a;^ + x'^z^. 

If we denote by Ui the subspace of Z^lx, y, z] generated by {51, 52, g?,,gi, 55}, then 
all the solutions with degree < 2 in each variable, are fi + Ui. 
Similarly we obtain: 

(1) /2 ~ X + and all the solutions are f2 + U2, where C/2 is the subspace 
generated by the polynomials 

hi^2 + x + xy + y'^, /12 = 2 + 2y + + 2y'^ , 

h3 = l + 2x + 2y^ + xy^, hi ^ y + 2y'^ + x^y, ^ 2y + y^ + x^y^. 

{2) = 1 + y + y"^ and all the solutions are /s + J/a, where C/3 is the subspace 
generated by the polynomials 

vi^2 + z + 2y + yz, V2 = l + 2z + 2y^ + y^z, 

V3 = 2 + z + 2y + yz'^, V4 = 1 + 2z + 2y'^ + y'^z'^, V5 = 2z + 2^. 
Finally one of the functions that can describe the genetic network is the following 
f{x, y, z) = {x + z + x'^ ,x + y"^ ,1 + y + y"^) 

4. Solution over the finite field GF(p") 

We can solve the problems DF(Zp) and (REP) using Lagrange interpolation over 
the field GF(p") = X 8, 9^. Let / : 5* Zp be a function with S" C Zp". Let 
151 = m be the cardinality of S. let {ai, a2, ■ ■ ■ , Q;„} be a fixed basis of K. There 
is a natural one to one correspondence between the sets Zp" and K , namely 

A : (ai, . . . , a„) 1-^ aiai H h a„a„. 

Let S — \{S) C K . Now we have the partially defined function / = A o / o A^^ : 
S ^%p. We denote the elements of S by a. 

Now, using the Lagrange interpolation formula we have the following: ai, . . ., 
are m distinct elements of the finite field K and /(ai) = 61, . . . , f{'a,n) — bm, 
with bi, . . . ,bm elements in Zp. We know that lip C K. We rewrite the problem 
DF(Zp) as follows: 

[DF(p"):] Let S C. K and let f : S K he a function. We want a polynomial 
P{x) e K[x] such that 

Pix) = f{x), for all ^ € ^ C is:. 

We can observe that this is the same problem (REP) if we consider S — {ri , . . . , r-m+i} 
and f{rj) = r^+i G K. So in the following we denote both problem by DF(p"). 

Using Lagrange Interpolation, we know that: there exists a polynomial P € K[x] 
of degree d < m — 1 such that P{ai) = e if for i = 1, . . . , m. The polynomial is 
given by 

m — 1 m — 1 

P{x)^^h W (a^ - afc)~^(a; - afe). 
Then, we have proved the following theorem. 
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Theorem 4.1. The problem DF{p"') has solutions over the field K = GF{p^^) using 
Lagrange Interpolation. That is, there exists a polynomial Pq G ^^[a;], such that 
P{){ai) = bi G K , for i = 1, . . . ,m. The degree of the polynomial Pq is less than or 
equal to m ^ 1. If I is the ideal of k[x] generated by P{x) = {x —'cii) ■ ■ ■ (x — Om); 
then all the solutions are given by Po{x) + G{x), where G{x) G /. 

Now, we have a new algorithm to solve the problem. We know by the Theorem 
14.11 that the solution is a polynomial of the following form: P{x) = J2T=o ^kx''- 
We evaluate P{x) in all the elements of S. Then we obtain a system of linear 
equations with one solution 

m — 1 

(//) Bko!^ = b, for i = 1, . . . , TO 

We want to remark that the system (//) has rank to, since ^ for i =^ j. 
Finally we have an output a polynomial in one variable with degree less than or 
equal to to — 1. 

Example 4.2. Let = {(0, 1), (1, 0), (1, 1), (2, 1)} C Zg^ and f : S ^ Z^. The 
function f has the following table of values. 



f{xi,X2) 





1 


2 





* 


1 


* 


1 


2 





* 


2 


* 


1 


* 



Let a be a root of the polynomial + X + 2 in Z3. Then — 2q! + 1, and 
a basis for G-F(3^) is {a, 1}. A natural correspondence is {xi,X2) ^ xia + X2 
We have ai = (0, 1), ai = 1; 02 = (1, 0), 02 = a; 03 = (1, 1), 03 — a + 1; and 
a4 = (2, 1), 04 = 2a + 1. Then, the particular solution is given by a polynomial 
Po{x) — cq + cix + C2X^ + c^x^. We evaluate in the four elements of GF(3^) using 
the table of values. The matrix of the system is the following: 

/II 1 1 1 \ 

1 a o? o? 2 

la + 12a + l 2a 

\12a + l 2 a + 2 1/ 

Solving the system we obtain the polynomial Pq{x) — + 2x + o^x^ + x^ . Now, we 
change the polynomial in two polynomials with two variables Xi, X2 and coefficients 
in Z3. We use the correspondence and obtain the following: 

f{xi,X2) = {2 + xi + 2xiX2 + x\,2 + 2xi + x^ + 2xiX2 + 2x\). 

5. Dynamical Systems over finite sets with different number of 

elements 

In this section we present a definition of dynamical system over finite sets. 

Definition 5.1. A dynamical system over finite sets is a pair (S*, /) such that 

(a) 5 C Zp". 

(b) A function, f = (/i, . . . , /„) : S ^ S, and each function fi : Si Zp is a 
partially defined function with Si C S for all i . 
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The DSF is a time discrete dynamical system, that is the dynamics is generated 
by iteration of the function /. 

Definition 5.2. The state space Sf of the dynamical system f : S S is a finite 
directed graph (digraph) with vertex set S and arc set A = {{x,y) S S : f{x) — y}, 
that is the ordered pair (x, y) stands for and arrow from a vertex x — . . . , Xn) S 
S to a vertex y = (t/i, . . . , ?/„) G S if and only if f{x) — y. 

Theorem 5.3. Let (S, f) be a dynamical system over finite sets. Then 

(a) f can he represented as a polynomial function in one variable over the finite 
field GF{p"), for some prime p. 

(b) If f = (/i: ■ ■ • J fn) then each function fi can he represented as a polynomial 
function in one variable or in several variables over Zp, for some prime p. 

(c) Part (a) and (b) hold for all prime number p > maXi\Si\, where \Si\ is the 
maximal number of different elements in the coordinate i of the set Si 

Proof. It is a consequence of the Sections 1 and 2. □ 



6. Examples and applications 

In this section we present two applications of the representation of the dynamical 
systems by polynomials over a finite fields. One very important things is to deter- 
mine the steady states, that is the elements x such that /(x) — x. We determine 
that in the first example. 

In the reverse engineering problem, we have a set of solutions and here we suggest 
a method for biologist to determine if one of the solution is the right one. 

In Fig. 1, an example of regulatory network is shown. This example of Gener- 
alized Logical Networks appear in ^5, • Here, we use the usual words for biologists. 
Gene 1 regulates genes 2 and 3, so that it has two thresholds (two values different 0) 
and the corresponding logical variable xi takes its value from {0, 1, 2}. Similarly, X2 
and X3 have one and two thresholds, respectively, and hence possible values {0, 1} 
and {0,1,2}. The functions /i(x), /2(x), and /3(x) need to be specified such as 
to be consistent with the threshold restrictions in the graph. Examples of logical 
functions allowed by the generalized logical method are shown in Fig. 1(b). 

Consider the case of /2(x). If xi ^ and ^3 ^ 0, so that xi and X3 have values 
above their first threshold, the inhibitory influences of genes 1 and 3 on gene 2 
become operative. Figure 1(b) indicates that X2 will tend to 0, that is, below the 
first threshold of the protein produced by gene 2. If either xi = or X3 = 0, that 
is, if only one of the inhibitory influences is operative, then gene 2 is moderately 
expressed. This is here represented by the value 1 for the image of X2. In general, 
several logical functions will be consistent with the threshold restrictions. Exactly 
which logical function is chosen may be motivated by biological considerations or 
may be a guess reflecting uncertainty about the structure of the system being 
studied. 

Example 6.1. 

G3 O 








2 


X2 





1 



Figure 1 (a) 
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f2{xi,X3) \ 1 2 \ I /3(xi,X3) \ 1 2 

111 2 2 1 

1 1 1 2 2 1 

2 \ 1 o\ I 2 \ 

Figure 1(h) 

Now, we describe this example using a dynamical system over finite sets. We 
have three genes, and the regulatory network is the following: 

3 O 

(1) The digraph: Y / j 

1 ^ 2 

(2) The variables x\ and 2:3 are in Z3 = {0, 1, 2}, and the variable is in Xi = 
{0,1}, 

(3) There arc several possibilities for functions /i and /2, but /s is the unique 
function. Then the polynomial functions /i, /2, and fs are: 

fl{x2) = 2X2 
f2{xi,X3) = 1 + 2x1x1 

fsixijXs) = 2 + xi + 2xz + x\X5 + 2x\ + x\ + 2x\xi 

+2xix1 + 2:^X3 

(4) The global function / : {Zs.f {Is.f, 

f(xi,X2,X3) = (fi(x2), ,f2{xi,X3), f3(xi,X3)). 

The state space has vertices V — I^^x X2 x Z3. We want to know the steady states 
of the dynamical system /, in general it is a very difficult problem. But, in this 
particular case we have three equations in three variables, 

2x2 = Xi, 1 + 2x\x\ = X2, 

2 + xi + 2x3 + X1X3 + 2x\ + x1 + 2x\x3 + 2xix\ + x\x\ = .T3 

For X2 we have only two values, if 2:2 = then x\ = Q and we obtain = 1 in 
the second equation, that is impossible. If a;2 = 1 then a;i = 2 and 1 + 2a;3 = 1 so 
2:3 = 0. We can check in the last equation and the only solution is (2, 1,0). 

Example 6.2. In the example that we present in Section 3, we have 'i^^ = 14, 348, 907 
different solutions. But, we can select the particular solution f{x,y, z) = {x + z + 
x^,x + y'^,l + y + y'^), and try to find which vectors in the state space go to the first 
state (1,2,0). 

Solving over Z3 the equations x + z + x"^ = 1, x + = 2, and 1 + y + j/^ = 0, 
we obtain that: 1 + y + y"^ = has one solution y = 1. So a; = 1 and z = 2 = 0( 
mod 2). Therefore, /(1, 1,0) = (1,2,0). Now, we can check that in the laboratory. 
Since the state space has only 18 elements, there are several functions with this 
property. 
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